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p£ : Abstract 

CNJ ' Current magnetic traps can be made so anisotropic that the atomic gas 



trapped inside behaves like quasi one or two dimensional system. Unlike the 
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• homogeneous case, quantum phase fluctuations do not destroy macroscopic 
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off-diagonal order of trapped Bose gases in d < 2. In the dilute limit, quan- 
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turn fluctuations increase, remain constant, and decrease with size for 3, 2, Id 



C$ • respectively. These behaviors are due to the combination of a finite gap and 
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i ■ the universal spectrum of the collective mode. 
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The technology of magnetic traps has advanced rapidly since the first series of reports 
of Bose-Einstein condensation in atomic gases of 87 Rb [|IJ, 23 Na 0, and 7 Li ||. The latest 
generation of magnetic traps ||] is capable of trapping 10 6 atoms in the ground state, 
three orders of magnitudes higher than that in the first Rb condensate experiment [pj. A 
common characteristic of these traps is that they are very anisotropic. The cylindrically 
symmetric harmonic potentials in these traps have transverse frequency u± about 20 times 
larger than the longitudinal frequency uj z . The resulting condensate looks like a cigar aligned 
along the symmetry axis z. With current technology, there appears no difficulty in achieving 
even higher asymmetry, say, with uj±/uo z ~ 100. In such limit, one can produce atomic gases 
with all the atoms lying in the lowest harmonic oscillator state in the xj/-plane, leaving the 



motion along z the only degrees of freedom. The system then behaves like a one dimensional 
(Id) Bose gas 0. Similarly, for uj±/uj z « 1, the system behaves like a 2d system [[7]. 

The possibility of realizing these low dimensional systems raises the question of the 
effects of quantum fluctuations. For homogeneous Bose systems, it is well known that phase 
fluctuations destroy finite (zero) temperature Bose condensation in dimension d < 2(< 1). 
The situation, however, is different for trapped Bose gases. For non-interacting bosons, 
finite spacing between the lowest and the next energy level allow BEC to occur at finite 
temperatures even in ID || |§ . The absence of gapless excitations also means that BEC will 
not be destroyed immediately as interactions between bosons are turned on. Nevertheless, for 
trapped bosons, quantum fluctuation effects due to interactions do show up in an unexpected 
way, which is the subject of this paper. We shall show that while the gapped nature of the 
collective modes eliminates the usual infrared divergence in all dimensions, the number of low 
energy modes contributing significantly to quantum fluctuation increases with the system 
size. As a result, quantum fluctuations suppresses but do not destroy the condensate. They 
also grow with size at a rate that increases with dimensionality. More precisely, we find that 

(I) The density matrix W(r, r') of a (^-dimensional dilute Bose gas in a harmonic trap 
exhibits macroscopic off-diagonal order of the form W(r,r')^ 0*(r)0(r') for |r — r'| of the 
order of the system size R, (both r and r' are inside the atom cloud); p a is the density of 
the system, 0(r) = y / p o (r)e _7d / 2 (up to a phase factor), and 7^ is a constant arising from the 
Gaussian fluctuations about the mean field solution : ji oc -R _1 ln(_R/cr), 7^ oc (R/a) d ~ 2 for 
d > 2, where a is the width of the ground state Gaussian of the d -dimensional harmonic 
well with frequency u, a = J%jMu3. These results show that within the dilute regime, the 
Gaussian approximation becomes more (less) accurate with increasing size in Id (3d). This 
anomalous demensionality dependence of the quantum fluctuations can be directly confirmed 
experimentally by studying how the momentum distribution depends on the number of 
particles. 

(II) The behavior in (I) is due to the gapped nature as well as the universal feature of the 
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collective mode. For d — 1, the collective modes are VL? = u y/£{£ + l)/2, I = 1,2,... For 



d = 2, they are il n , m = ujy(n 2 — m 2 + 2n)/2, n = \m\ + 2p, p = 0, 1,2, .., m =integer. 

We shall first derive these results, and then discuss their implications for real atomic 
gases. Let us consider a ci-dimensional dilute Bose gas with a partition function Z = 
Tre _/3(ff_/iAr) , where /x is the chemical potential, K = H - pN= J d d v[(h 2 /2M) • 
Vijj+(U - p) +(g d /2)-tp + -ijj + -tjjtjj], g d > 0, U(r) = (l/2)Muj 2 \r\ 2 is a d -dimensional 

harmonic well, and V is a rf-dimensional gradient. H is an effective low energy Hamiltonian 
where the Fourier component of tp with wavelength shorter than an atomic scale a c has been 
coarse-grained out. 

It is convenient to write the partition function as a coherent state path integral Z = 

J D(j)*D(j)e~ s l h , where 
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WW + ^T|V0| 2 + (U - /i)|0| 2 + ^|0| 4 



(1) 



2M 1 ri v 2 
with / = / d^r Jo' 3 dr. Anticipating a condensate, we write <p = e t8 y/p . The stationary point 
is a constant phase O and a density p (r) satisfying the Gross-Piteavksii (GP) equation, 

h 2 



2M 



V sJVo + (U- fj)y/fa + gdpo^/Vo = o. 



(2) 



Expanding S about (9 , p (r)) to quadratic order in phase fluctuations Q\ and density 
fluctuations p 1 , (pi — p — p Q , Q\ = 6 — o ), we have S = S [p , 9 Q ] + Si, where 

Si = J \ih P id T 6i + ie 1 -f-e 1 + ^pi-d- P i~\ (3) 

e 1 -f-e l = -6iV ■ (p [h 2 /M] V0i) (4) 
Pi-G-pi — g d p\ - (h 2 /2M)[( Pl /^p- )V 2 ( Pl /^p- ) - (pVpTV 2 ^)] (5) 

One can also integrate out the density fluctuations and obtain an effective action of phase 
fluctuations 6>i, 



Vase = / [^ 2 (d T ei) ■ g- 1 ■ {d T ei) + -e 1 .f-e 1 



(6) 



Quantum fluctuations are caused by zero point motions of the excitations. The low- 
lying excitations are phonons, corresponding to density oscillations inside the atomic cloud. 



At high excitation energy, the period of density oscillations become comparable to cloud 
size, the excitations leak out of the cloud and become single-particle like, and cannot be 
adequately described by the phase action eq. The equations of motion of the collective 
modes in real time t can be obtained by first analytically continuing the action eq.(^) to 
real time (r — > it) and then finding the stationary conditions. They are Tid t 6\ = Gpi and 
—hdtpi = T9\. The wavefunctions and the frequencies of the collective modes (labelled by 
a set of quantum numbers a = {««}) are the eigenfunctions and eigenvalues of the operator 
GT, i.e. 



GTu a = h 2 tt 2 a u a . (7) 

The Q = mode (corresponding to a constant phase) is excluded from the collective excita- 
tion spectrum because it is simply the ground state. Since all other modes are orthogonal 

to the Q = mode, they have zero spatial average a fact that shall be of use to us later. 

In 1956, Penrose and Onsager || pointed out a unified description of Bose-Einstein 
condensation (BEC) for both noninteracting and interacting systems based on the general 
structure of the density matrix ^(r, r') =< tp + (r)tp(r f ) >= I] a A Q ,z/*(r)i/ Q (r') , where X a 
and u a are the eigenvalues X a and normalized eigenfunctions of W. BEC is characterized |5| 
by the fact that below a transition temperature T c , one of the eigenvalues (say Ao) becomes 
order N compared with all others, 

W = X u*(r)u (r') + K(r,r'), A /A^ ~ N, (8) 

where K denotes the sum of those eigenstates A^o an d (for homogeneous systems) has the 
property that K(r,r') — > as |r — r'| — > oo. The first term of eq.(|8]), which represents 
condensation, is referred to off-diagonal- long-range order |H|. For trapped Bose gases, 



however, both r and r' must remain inside the atomic cloud. The large distance behavior 
of W is masked by finite size effects. In spite of this, the appearance of a macroscopic 
eigenvalue as described in eq. (|8[) remains a proper characterization of BEC || . 

To demonstrate the existence of macroscopic off-diagonal order in trapped Bose gases in 
all dimensions at T = 0, we write 



W(r, r') =< v /p(r)p(r')e- i ^W- ( r ')] >= ^ Po {v) Po {v') [er F ^ + F(r, r')) , (9) 
where F(r,r') = \ < [9\{v) - #i(r')] 2 >, and 



2 \\Po(r) Mr')/ , 



(pi(r)pi(rQ) _ <pf(r)> _ (p 2 (r')) + Q ( ( Pi\ 
4p (r)p (r') 8p 2 (r) 8p 2 (r') Up oy J 



3/2N 



(10) 



As we shall see later, the term Y(r, r') is smaller that the first term in eq.@ by at least a 
factor of (cx/R) 4 where R is the size of the cloud. We can then ignore the term K(r, r') in 
eq.(||). Using eq.(^), it is straightforward to show that 

F(r, r') = Y, \< r\G 1/2 \a >-< r'lG^a > 2 — (11) 

where < r\at >= u a (r). At first sight, the phase fluctuation appears to decrease with the 
size of the system R because w 2 scales as R~ d . This, however, is not true because the ot- 
sum (which depends on the number of excitations contributing significantly to the phase 
fluctuation) is also R dependent. To evaluate F, we need to study the collective modes in 
greater detail. For systems with large number of particles, it is known that Thomas-Fermi 
approximation (TFA) becomes accurate | IT| . TFA ignores the gradient term in eq. (0). It 



implies that 

p (r) = ^^(l-^6(i? 2 -r 2 ), p^-M^R\ (12) 

where Q(x) = 1 (or 0) of x > (or < 0). R is the size of the system. The number constraint 
N = J p implies that the system size R is related to iV as 

R/a= [g d N(d + 2)/(huJC d cr d )] 1/{d+2) , (13) 

where Qcr d is the volume of a rf-dimensional sphere of radius a. The fact the p Q oc 
gj 1 in TFA means that T oc g d l . In addition, TFA also implies that the term 
(fr 2 /4M)p 2 /p^/ 2 V 2 v /p^)in G (eq.(^j)) is negligible compared to g&. Furthermore, for density 
oscillations with period small compared to cloud size (phonon regime, see below), (h 2 /AM) 



[(pi/ v /p^)V 2 (pi/ y/fo)] << Qdi an d can also be ignored. As a result, G = ga, and the explicit 
g dependence in eq.(0) is cancelled out. For harmonic potentials, eq.(^) becomes Ul2f 

^ 2 V([i? 2 -r 2 ]Vn Q ) =n 2 a u a , (14) 

For d — 1, eq. (|i4|) is the Legendre equation. The normalized eigenfunctions are ug(r) = 
(2R)~ 1 / 2 Pe(r / R), where Pe(x) are Legendre polynomials. The frequency of is = 
(a?/ \fi)\jt{t + 1), £=1,2,... In d = 2, it is straightforward to show that the eigenfunctions 
are of the form w n)m (r,0) = e~ im ^ Y7e=\ m \ Ce(r / R) e with n = \m\ + 2p,p = 0,1,2,... The 
corresponding frequencies are VL nrn = (&/ \f2)\Jn 2 — m 2 + 2n. We have thus established 



statement (II). The collective mode in d = 3 has been derived by Stringari fT2]| . Because 
of the gapped nature of the collective mode, the phase fluctuation eq.(|Tl~D reduces its zero 
temperature form for ksT < Two. In other words, the existence of macroscopic off-diagonal 
order at T = immediately implies its existence at sufficiently low temperatures. We shall 
therefore from now on focus on the T = case. 



Dilute condition and phonon cutoff : To evaluate eq. (|Tl|) , it is necessary to discuss a 
number of implicit constraints implied by the low energy model (eq. ([[])), as well as how the 
collective modes behave as their energy increases. First of all, to stay in the dilute limit, the 
density at the center of the trap must be less that the inverse volume of the atomic cutoff, 
p o (0) < (crfaf) -1 . This limits the size of the cloud to be smaller than R m , 

R\^ gd f Rm\^ 

) < Tt h = ( ~~ ) ■ dilute condition. (15) 



For a cloud of radius R satisfying eq.([L5|), the collective modes will be labelled by d quan- 
tum numbers ctj, i = l,..d, denoting the number of nodes in d different curvilinear di- 
rections. The spatial variations of this mode along a curvilinear direction is given by a 
averaged wavevector ki ~ oti/R. The corresponding gradient energy is h 2 kf/2M. As the 
energy of the collective mode increases (i.e. increasing «j), ki will increase to the point 
where the condition (h 2 /4M) [(/Oi/y / p^)V 2 ( / oi/ v /p^)] << gapf breakdowns. This occurs 
when h 2 k 2 /2M ~ g d p(0) = \Muj 2 R 2 , which is h = R/a 2 , or 

6 



<y.i = (R/cr) = a p h : Phonon cutoff (16) 

Obviously, both R m and a p h are only defined up to some numerical factor. The collective 
modes with quantum numbers below a p h are described by eg. . For frequencies Q > 
£la ph , the excitation eq.(|7] ) will not be phonon like. Instead, they leak out of the cloud 
and become single-particle like. The action eq.([5]) is no longer sufficient. Eventually, the 
average wavevector of the collective mode will approach the "high energy" cutoff ki ~ l/a c , 
corresponding to 

a.i ~ R/a sc = a*, : high energy cutoff, (17) 

in which case the low energy model eq. ([!]) is no longer valid. Excitations between the phonon 
cutoff a p h and the high energy cutoff a* will be referred to as "high energy" modes. (Clearly 
we need to have R/a c > (R/cr) 2 , (i.e. cr/a c > -R/cr) for the cutoffs in eq.fll6|) and eq.flT^) to 
be consistent. For a/a c ~ 10 2 , this condition is easily satisfied even for large clouds with 
R/cr ~ 20. 

Now we return to eq. @. We note that the first term in eq.flTOD is odd in 9 and vanishes 
by symmetry. For the second term in eq.fllQl), it is straightforward to show (13| that its 



contribution from a given phonon mode of average momentum k ~ 1/R is down from the 
corresponding phase fluctuation by h 2 k 2 / (Mg d p ) ~ (cr/R) 4 1 for large clouds. We 
can therefore ignor the density fluctuation term etc in eq.(|9D as mentioned before. Writing 
F(r,r')= +J2( R /a)i <ai <R/a c )(•••) = F P h + Rhigh, where (...) is the summand of 

eq. (|TTD . The leading R dependence comes from F ph . The physics of F high and high energy 
density fluctuations is analogous to the similar high energy contributions in homogeneous 
systems, which are known to be less important than the phonon contributions for large 
|r — r'| ||13|| . By noting that within the phonon regime, (i.e. for excitations with «j < a p h), 
the operator G reduces to the constant ga, Eq.flTTD implies that 



F»M= E ^-K(r)-«a(r')| 2 . (18) 
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Denoting the spatial average and fluctuation of a function /(r) as / and A(/), i.e 
/ = / + A(/), we have 

1 



F ph (r, r') = 7d + r? d (r, r'), 7d 



Cd -^ d a<(i?/cr) 2 ^ a 



E 



(19) 



and %(r,r') = (g d /2hn a ) £ a<W(T)a [A(K(r)| 2 + |^(r')| 2 ) -(u a (r)*u a (r') + c.c.)]. 

To estimate we note that there is a one-to-one correspondence between the solutions 
of eq. ([14]) and those of the equation 



u 2 R 2 /2)V 2 u = tt 2 u, 



(20) 



as both of them are labelled by the same quantum numbers. Moreover, the solutions of 
eq.flHD resemble those of eq. Q20|) as the number of nodes increases so that the gradient 
term becomes more important. We thus approximate collective mode frequencies and wave- 
functions in eq. fll8l) by those of eq. (p0|) . The normalized functions u a (r) will be replaced 
normalized plane waves Uk(r) = e lk r / \/c d R d , where k is a d dimensional wave- vector. kiR is 
the correspondence of on as it counts the number of nodes along the z-th (linear) direction. 
The frequency f2 Q in eq . ([T8|) corresponds to cu(Rk) j\f2 in eq . (|20 ) . With this approximation, 
eq.(|l8D becomes 

1 y/2g d [W°) 2 d d (kR) ( 1 - cosk • (r - r') \ 



F ph (r,r') 



(21) 



c d R d hoo h (27r) d V kR 

The constant and the cosine term in the bracket in eq.(pl|) correspond to the average and 
fluctuation term (7^ and rj d in eq.(|19|) respectively). The lower limit is due to the fact that 
the lowest collective mode has an average wavelength R (or wavevector 1/R). 
Eq.(J2Tp implies that 

R^ 



yl ( 91 ) (- 

7T \2ahuj) \R 



21n 



a 



Ci ( ^] + Ci(2x) 



(22) 



4vt KhuJTia 2 J \R 



a 



R^ 2 
a . 



r(l/2) [J (a)H_ 1 (a)-H(a)J_ 1 (a)]H 



a=xR 2 /(7 2 



(23) 
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where x = |r — r'\/R, Ci is the Cosine integral, J a , J_i are Bessel functions, T is the Gamma 
function, and H Q , H_i are the Struve functions. The first and second terms in these expres- 
sions are the constant 7^ and the function rjd in eq. (|T5D . It is easy to show that rjd — > 
as |r — r'| — > 00. Even for |r — r'\/a ~ 2, rjd << jd provided R/a > 5. In other words, 
for large clouds, we have F ph ~ 7^ when |r — r'| ~ R/2. Eq.@ and (|TT|), and the fact 
that F ~ F p h implies statement (I). Note that although the cutoff a p h in our derivation 
contains an uncertain numerical factor, it does not affect the R dependence of our results. 
Our results also leads to the interesting conceptual point that for a true Id system with a 
5 -function potential of finite strength, the mean field off-diagonal form of density matrix 
becomes exact as R — > 00. 

The anomalous R dependence of our results can be qualitatively understood simply as 
follows. The density of trapped bosons is not arbitrary, but increases with size R according 
to the scaling relation p ~ R 2 . For a bulk homogenous system, the dimensionless coupling 
constant is ^f^p( d ~ 2 )/ 2 ; and so increasing density implies stronger (weaker) coupling in 3D 
{ID). Indeed, the R dependence of the condensate depletion in d— dimension is obtained 
if we replace the trapped boson system with a homogenous system of size R and density 
p ~ R 2 . 

Atomic gases in quasi Id and 2d limit : Our previous discussions show that atomic 
Bose condensates exist in these quasi low dimensional regimes at sufficiently low (non- 
zero) temperatures. These condensates can be detected through their collective modes in 
Statement (II). In addition, the momentum distribution n(k) (i.e. Fourier transform of W) 
will reflect the size dependence of the quantum fluctuation. 

When only the lowest harmonic quantum state in the tightly confining direction is oc- 
cupied, the coupling constant in the effective quasi Id and 2d theory is given by 

gi = 2hu_ L a sc , g 2 = 2\/27ihuj z a z a sc , (25) 

where a z and a± are the widths of the ground state harmonic oscillator along z and in the 

9 



xy-plane, and a sc is still the s-wave scattering length in three dimension. The quantum 
suppression parameters 7^ read from first terms of Fp h (r, r'), d — 1, 2, 3, are then 

(26 > 

From eq. fll3"D , we see that 7^ ~ (a sc ) ( '' i+1 ^ ( -' i+2 ' ) , and so these effects are particularly notice- 
able in systems with large a sc . 

In summary, we have shown that the quantum fluctuations in trapped boson systems ex- 
hibit anomalous dimensionality dependence. Our results may be observeable in experiments 
on quasi-l-D and 2D systems with anisotropic traps. 
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